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Optical properties of black hole in the presence of plasma: shadow 
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We have studied photon motion around axially symmetric rotating Kerr black hole in the presence 
of plasma with radial power-law density. It is shown that in the presence of plasma the observed 
shape and size of shadow changes depending on i) plasma parameters, ii) black hole spin and iii) 
inclination angle between observer plane and axis of rotation of black hole. In order to extract pure 
effect of plasma influence on black hole image the particular case of the Schwarzschild black hole 
has also been investigated and it has been shown that i) the photon sphere around the spherical 
symmetric black hole is left unchanged under the plasma influence, ii) however the Schwarzschild 
black hole shadow size in plasma is reduced due to the refraction of the electromagnetic radiation 
in plasma environment of black hole. The study of the energy emission from the black hole in 
plasma shows that in the presence of plasma the maximal energy emission rate from the black hole 
decreases. 

PACS numbers: 04.20.-q; 95.30.Sf; 98.62.Sb 


I. INTRODUCTION 

The study of astrophysical processes in plasma medium 
surrounding black hole becomes very interesting and im¬ 
portant due to the evidence for the presence of black 
holes at the centres of the galaxies [1-3]. For exam¬ 
ple, the gravitational lensing in inhomogeneous and ho¬ 
mogeneous plasma around black holes has been recently 
studied in [ - ] as extension of vacuum studies (see, e.g. 
[9, 10]). 

From the literature it is known that the black hole 
shadow is appeared by the gravitational lensing effect, 
see, e.g. [ 1 , 11-13]. If black hole is placed between a 
bright source and far observer, dark zone is created in 
the source image by photon fall inside black hole which 
is commonly called shadow of black hole. Recently, this 
effect is investigated by many authors for the different 
black holes (see, e.g. [ - ]). The silhouette shape of 

an extremely rotating black hole has been investigated 
by Bardeen [18]. Our previous studies on the shadow 
of black hole are related to the non-Kerr [ I], Horava- 
Lifshitz [T ], Kerr-Taub-NUT [If] and Myers-Perry [20] 
black holes. A new coordinate-independent formalism for 
characterization of a black-hole shadow has been recently 
developed in [21]. 

Shape of black hole is determined through boundary 
of the shadow which can be studied by application of the 
null geodesic equations. The presence of plasma in the 
vicinity of black holes changes the equations of motion 
of photons which may lead to the modification of black 
hole shadow by the influence of plasma. In this paper 
our main goal is to consider silhouette of shadow of axi¬ 


ally symmetric black hole using the equations of motion 
for photons in plasma with radial power-law density. We 
would like to underline that very recently, influence of a 
non-magnetized cold plasma with the radially dependent 
density to black hole shadow has been studied in [ ] us¬ 
ing the different alternate approach. In addition Rogers 
[ ] has studied the photon motion around black hole 

surrounded by plasma. 

The paper is arranged as following. In Sect. II, we 
consider the equations of motion of photons around axi¬ 
ally symmetric black hole in the presence of plasma. In 
Sect. Ill we study the shadow of the axial-symmetric 
black hole in the presence of plasma. As particular case 
in subsections III A and IIIB we study the shadow and 
the energy emission from the sphericaly symmetric black 
hole. Finally, in Sect. IV we briefly summarize the results 
found. 

Throughout the paper, we use a system of geometric 
units in which G = 1 = c. Greek indices run from 0 to 3. 


II. PHOTON MOTION AROUND THE BLACK 
HOLE IN THE PRESENCE OF PLASMA 


The rotating black hole is described by the space-time 
metric, which in the standard Boyer-Lindquist coordi¬ 
nates, can be written in the form 
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where as usual M and a are the total mass and the spin 
parameter of the black hole. 

In this paper we will consider the plasma surrounding 
the central axially symmetric black hole. The refraction 
index of the plasma will be n = n(x l ,uj) where the photon 
frequency measured by observer with velocity u a is u). In 
this case the effective energy of photon has the form fuj = 
—p a u a . The refraction index of the plasma as a function 
of the photon four-momentum has been obtained in [ ] 

and has the following form: 


n 2 = 1 + 


PaP a 

(ppu p ) 2 ’ 
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and for the vacuum case one has the relation n = 1. The 
Hamiltonian for the photon around an arbitrary black 
hole surrounded by plasma has the following form 


changes with the radial position and is redshifted in the 
spacetime. Assume the Killing vector to have compo¬ 
nents 

r = (m,0,o[) ; & = </oo( —1,0,0, o) , (8) 


so that = fco =const. The frequency of an electromag¬ 
netic wave emitted at radial position r and measured by 
an observer with four-velocity u a {l/y/—goo,0, 0, 0} par¬ 
allel to (i.e. a static observer) will be governed by the 
following equation 


V~ #oo^(t) = = const . 


(9) 


One may introduce a specific form for the plasma fre¬ 
quency for analytic processing, assuming that the refrac¬ 
tive index has the general form 

n 2 = l-4> (10) 

where uj e is usually called plasma frequency. Now using 
the Hamilton-Jacobi equation which defines the equation 
of motion of the photons for a given space-time geometry 
[4, 22, 24]: 
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where p a = dS/dx a . Using method of separation of vari¬ 
ables the Jacobi action S can be written as [ 3, 23]: 
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Following to the derivation of a gravitational redshift 
discussed in [25] we will assume that the spacetime sta¬ 
tionary allows existence of a timelike Killing vector 
obeying to the Killing equations 


Car,/3 T C/3;ck — 0 • (5) 

Then one can introduce two frequencies of electromag¬ 
netic waves using null wave-vector k a the first one is the 
frequency measured by an observer with four-velocity u a 
and defined as 


uj = —k a u a , (6) 

while the second one is the frequency associated with the 
timelike Killing vector and defined as 


where £, £ are conservative quantities as angular mo¬ 
mentum and energy of the test particles. 

For trajectories of the photons we have the following 
set of the equations: 

E-^- = a(C — n 2 £a sin 2 6) 
da 
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can be derived from the Hamilton-Jacobi equation, where 
the functions 7 Z(r) and ©(#) are introduced as 


^ = -k a £ a . (7) 

The frequency (6) depends on the observer chosen and 
is therefore a function of position, while the frequency (7) 
is a conserved quantity that remains unchanged along the 
trajectory followed by the electromagnetic wave. One 
can apply this property to measure how the frequency 


n = 

0 = 


[(r 2 + a 2 )£ — aC\ 2 + (r 2 + a 2 ) 2 (n 2 — 1)£ 2 
-A [K + (£ — a£) 2 ] , (17) 

JC + cos 2 6 (a 2 £ 2 -~ 2 — 

\ sim 0 

— (n 2 — l)a 2 S 2 sin 2 6 , 


( 18 ) 


















3 


and the Carter constant as JC. 

For calculation examples one needs the analytical ex¬ 
pression of the llama frequency uj e which for the electron 
plasma has the following form 
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where e and m e are the electron charge and mass re¬ 
spectively, and N(r) is the the plasma number density. 
Following the work by Rogers [ 22 ] here we consider a 
radial power-law density 


«(r) = ( 20 ) 

where h > 0 , such that 
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As an example here we get the value for power h as 1 [ 2 < ]. 
For this value we plot the radial dependence of the effec¬ 
tive potential V e & of the radial motion of the photons 
defined as 
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The radial dependence of the effective potential for dif¬ 
ferent values of plasma refraction n and black hole spin 
a has been presented in Fig. 1. In the Fig.l the left 
plot corresponds to the case when refraction parameter 
of the plasma is n 2 = 0.2 , 0.44 , 0.89 (dotted, dashed 
and solid lines, respectively) at the position r = 3M; 
middle plot corresponds to the case when the refraction 
parameter is n 2 = 0.19 , 0.42 , 0.88 corresponding to 
dotted, dashed and solid lines, respectively, at the posi¬ 
tion r = 3M; right plot represents the radial dependence 
of the effective potential when refraction parameter is 
n 2 = 0.14 , 0.39 , 0.88 corresponding to dotted, dashed 
and solid lines, respectively, at the position r = 3 M. 


III. SHADOW OF BLACK HOLE IN THE 
PRESENCE OF THE PLASMA 

In this section we consider the shadow cast by black 
hole surrounded by plasma. If black hole surrounded by 
plasma originated between the light source and the ob¬ 
server, then the latter can observe the black spot on the 
bright background. The observer at the infinity can only 
observe the light beam scattered away and due to cap¬ 
turing of the photons by the black hole the shaded area 
on the sky would be appeared. This spot corresponds to 
the shadow of the black hole and its boundary can be 
defined using the equation of motion of photons given by 
expressions (13)-(16) around black hole surrounded by 
plasma. 

In order to describe the apparent shape of the the black 
hole surrounded by plasma we need to consider the closed 


orbits around it. Since the equations of motion depend 
on conserved quantities £, C and the Carter constant /C, 
it is convenient to parametrize them using the normalised 
parameters £ = £/£ and 77 = /C/£ 2 . The silhouette of 
the black hole shadow in the presence of the plasma can 
be found using the conditions 

lZ{r) = 0 = dlZ(r)/dr. 

Using these equations one can easily find the expressions 
for the parameters £ and 77 in the form 
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where we have used the following notations 
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and prime denotes the differentiation with respect to ra¬ 
dial coordinate r. 

The boundary of the black hole’s shadow can be fully 
determined through the expessions (23)-(24). However, 
the shadow will be observed at ’observer’s sky’, which 
can be referenced by the celestial coordinates related to 
the real astronomical measurements. The celestial coor¬ 
dinates are defined as 
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a — lim —r 0 sm 0 0 — 

r 0 —5>oo Y dr 

o y 2 d0 

P = lim r 0 — . 

r 0 —>-oo dr 


(28) 

(29) 


Using the equations of motion (13)-(16) one can easily 
find the relations for the celestial coordinates in the form 
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(31) 


for the case when black hole is surrounded by plasma. 

In Fig 2 the shadow of the rotating black hole for the 
different values of black hole rotation parameter a, in¬ 
clination angle #o between the observer and the axis of 
the rotation is represented. In this figures we choose the 
plasma frequency in the form uj e /uj% = k/r. From the 
Fig. 2 one can observe the change of the size and shape of 
the rotating black hole surrounded by plasma. Physical 
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FIG. 1: The radial dependence of the effective potential of radial motion of photons for the different values of rotation 
parameter and refraction index of the plasma. Here the quantity V e s is normalised by the energy of the photon S. 


reason for this is due to gravitational redshift of photons 
in the gravitational field of the black hole. The frequency 
change due to gravitational redshift affects on the plasma 
refraction index. 


A. Shadow of non-rotating black hole 


Now in order to extract pure plasma effects we will 
concentrate at the special case when the black hole is 
non-rotating and the size of the black hole shadow can 
be observed (see, e.g. [ ]). In the case of the static black 
hole the shape of the black hole is circle and the radius of 
the shadow will be changed by plasma effects. Using the 
expressions (30) and (31) one can easily find the radius 
of the shadow of static black hole surrounded by plasma 
in the form: 


Rsh — 


2 r 6 {r — M)n z + rrm (r — M) 


n(r — M ) 

—2 r 2 M 2 + 2Mr 2 |nr 2 (n + rn f ) — (4 n + 3 rn') 
xnMr + M 2 (l + 3 n 2 + 2rrm / )| 


, (32) 


here r is the unstable circular orbits of photons de¬ 
fined by dr/da = 0 and dV e ^/dr = 0. In the absence 
of the plasma one has the standard value of the pho¬ 
ton sphere radius as r = 3 M and for shadow radius 
R s h = 3^3 [ 26 , 27]. In the presence of the plasma we 
will have different value for the photon sphere radius and 
consequently different shadow radius of the boundary of 
black hole shadow. In the Fig. 3 the dependence of the 
radius of shadow of the static black hole from the plasma 
parameters has been presented which shows that the ra¬ 
dius of the shadow of black hole surrounded by inhomo¬ 
geneous plasma decreases. It is similar to the results of 
the paper [7]. 


B. Emission energy of black holes in plasma 


For the completeness of our study here we evaluate 
rate of the energy emission from the black hole in plasma 
using the expression for the Hawking radiation at the 
frequency ft as [28, 29] 

d?E(Q) _ e lr 2 aii rn 3 
dttdt “expU/T-1 5 [ J 


where T = k/2i\ is the Hawking temperature, and k is 
the surface gravity. Here, for simplicity we consider the 
special case when the black hole is non rotating, and the 
background spacetime is spherically-symmetric. 

At the horizon the temperature T of the black hole is 


T = 



(34) 


The limiting constant (Tu m 

Olim ~ ^Rsh 


defines the value of the absorption cross section vibration 
for spherically symmetric black hole and R s h is given by 
expression (32). 

Consequently, one can get 

*E(Sl) 27T 3 i4 3 

dtddt e n / T — 1 

that the energy of radiation of black hole in plasma de¬ 
pends on the size of its shadow. 

The dependence of energy emission rate from fre¬ 
quency for the different values of plasma parameters uo e 
is shown in Fig. 4. One can see that with the increasing 
plasma parameter uo e the maximum value of energy emis¬ 
sion rate decreases, caused by radius of shadow decrease. 


IV. CONCLUSIONS 

In this paper we have studied shadow and emission 
rate of axial symmetric black hole in presence of plasma 
with radial power-law density. The obtained results can 
be summarized as follows. 
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FIG. 2: The shadow of the black hole surrounded by plasma for the different values of the rotation parameter a, inclination 
ante between observer and the axis of the rotation #o, and the refraction index. The solid lines in the plots correspond to the 
vacuum case, while for dashed lines we choose the plasma frequency uJe/w^ — k/r and ( k/M ) 2 = 0.5 . 


• In the presence of plasma the observed shape and 


size of shadow changes depending on i) plasma pa- 
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FIG. 3: The dependence of the Schwarzschild black hole 
shadow on the plasma frequentcey parameter. Here we take 
the plasma frequency as co e /uj^ — k/r. 



FIG. 4: Energy emission from black hole for the different 
values of k/M : Solid line corresponds to the vacuum case 
(k/M = 0), Dashed line corresponds to the case when k/M = 
0.4, dot-dashed line corresponds to the case when k/M = 
0.6, dotted line corresponds to the case when k/M = 0.8. 
Here we take the plasma frequency as co e /uj^ — k/r and £7 is 
normalised to the T. 

rameters, ii) black hole spin and iii) inclination an¬ 
gle between observer plane and axis of rotation of 


black hole. 

• In order to extract pure effect of plasma influ¬ 
ence on black hole image the particular case of the 
Schwarzschild black hole has also been investigated. 
It is shown that under influence of plasma the ob¬ 
served size of shadow of the spherical symmetric 
black hole becomes smaller than that in the vac¬ 
uum case. So it has been shown that i) the photon 
sphere around the spherical symmetric black hole 
is left unchanged under the plasma influence, ii) 
however the Schwarzschild black hole shadow size 
in plasma is reduced due to the refraction of the 
electromagnetic radiation in plasma environment of 
black hole. 

• The study of the energy emission from the black 
hole in plasma has shown that with the increase of 
the dimensionless plasma parameter the maximum 
value of energy emission rate from the black hole 
decreases due to the decrease of the size of black 
hole shadow. 

In the future work we plan to study shadow and re¬ 
lated optical properties of different types of gravitational 
compact objects in the presence of plasma in more detail 
and in more astrophysically relevant cases. 
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